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ABSTRACT 

The  Oriented  Scintillation  Spectrometer  Experiment  (OSSE)  on  board  the 
Compton  Gamma  Ray  Observatory  detected  the  Vela  Pulsar  (PSR  B0833-45)  during 
August- September  1991,  April-May  1992,  and  August  1993.  Observed  lightcurves 
have  a  two-peak  pulse  profile  similar  to  that  observed  at  higher  energies,  although 
the  second  peak  may  be  wider  in  the  OSSE  lightcurve.  Pulsed  emission  in  the  first 
gamma-ray  peak  was  detected  with  4.6<r  statistical  significance  in  the  0.07-0.6  MeV 
band  in  the  sum  of  all  three  observing  periods.  The  second  gamma  ray  peak  was 
detected  at  no  more  than  3<r  significance  in  the  same  band.  Due  to  the  low  statistical 
significance  of  the  observations,  little  can  be  said  concerning  longer-term  temporal 
variability.  The  spectrum  is  hard  at  lower  energies  and,  in  combination  with  higher 
energy  data,  appears  to  require  a  break  in  the  20  MeV  region.  OSSE  also  observed 
Geminga  during  July  1992,  December  1993,  and  July  1994.  No  significant  pulsed  or 
time  averaged  emission  was  observed  on  any  occasion.  Upper  limits  to  the  pulsed 
emission  suggest,  but  do  not  require,  a  break  from  the  extrapolation  of  the  spectrum 
measured  at  higher  energies. 


Subject  headings:  gamma  rays:  observations  -  stars:  pulsars:  individual  (Vela  Pulsar, 
Geminga) 


1.  Introduction 


The  number  of  known  gamma-ray  pulsars  has,  over  the  past  few  years,  increased  from  two 
(Crab  and  Vela)  to  seven  with  the  addition  of  PSR  B1509-58,  PSR  B1706-44,  PSR  B1055-52,  PSR 
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B1951+32  and  the  discovery  that  the  gamma-ray  object  Geminga  is  a  pulsar  (Halpern  &  Holt, 
1992).  The  Oriented  Scintillation  Spectrometer  Experiment  (OSSE)  on  board  Compton  Gamma 
Ray  Observatory  (CGRO)  has  previously  reported  detections  of  the  Crab  Pulsar  (Ulmer  et  al. 
1994)  and  PSR  B1509-58  (Matz  et  al.  1994).  We  report  here  on  OSSE  observations  of  Vela  and 
Geminga.  Neither  of  these  objects  is  known  as  a  bright  hard  X-ray  or  low  energy  gamma-ray 
emitter,  although  both  are  relatively  nearby  and  both  are  intense  sources  of  gamma  rays  above 
100  MeV. 

The  Vela  Pulsar  (PSR  B0833-45),  with  a  period  of  0.089  seconds,  is  the  brightest  extrasolar 
object  in  the  gamma-ray  sky  at  energies  above  100  MeV.  It  has  been  observed  extensively  at  these 
energies  by  SAS  II  (Thompson  et  al.  1977),  COS-B  (Bennett  et  al.  1977;  Grenier,  Hermsen  & 
Clear  1988)  and  EGRET  (Kanbach  et  al.  1994).  In  addition  to  radio  emission,  the  Vela  Pulsar 
also  emits  faint  optical  (Wallace  et  al.  1977)  and  soft  X-ray  pulsations  (Ogelman,  Finley  & 
Zimmermann  1993).  While  the  radio  lightcurve  exhibits  a  single  peak,  the  gamma-ray  lightcurve 
is  double-peaked,  with  the  first  peak  lagging  the  radio  peak  by  ~  0.11  in  phase  and  the  two 
gamma-ray  peaks  separated  in  phase  by  ~  0.42. 

The  pulsed  X-ray  spectrum,  detected  by  ROSAT,  was  initially  reported  (Ogelman  et  al.  1993) 
to  be  consistent  with  thermal  emission  at  or  near  the  neutron  star  surface  rather  than  emission 
from  an  cascade  in  the  magnetosphere,  the  presumed  emission  mechanism  at  other  observed 
frequencies.  A  more  recent  interpretation  (Ogelman  1993)  suggests  that  the  total  spectrum  from 
the  point  source  co-located  with  the  pulsar  can  also  be  modeled  by  soft  and  hard  components, 
the  latter  of  which  may  indicate  a  magnetospheric  contribution.  However,  the  pulsed  spectrum  is 
still  best  represented  by  a  soft  blackbody  spectrum  and  is  primarily  observed  below  1  keV.  The 
lightcurve  at  soft  X-ray  energies  consists  of  a  broad  feature  (suggested  by  Ogelman  (1993)  to  be 
a  double  pulse)  lagging  the  radio  pulse  by  approximately  0.5  in  phase.  The  centroid  of  the  broad 
X-ray  peak  leads  the  second  gamma  ray  peak  as  shown  by  Kanbach  et  al.  (1994)  by  ~  0.1  in 
phase  (see  Figure  2). 

Until  recently,  Vela  has  been  an  elusive  target  in  the  hard  X-ray  and  low-energy  gamma-ray 
bands.  An  observation  of  0.3-30  MeV  pulsed  emission  was  reported  by  Turner  et  al.  (1984)  31 
days  after  a  major  pulsar  period  glitch,  but  negative  results  have  been  reported  by  Ulmer  et  al. 
(1991)  using  data  from  the  gamma-ray  spectrometer  aboard  HEAO  3,  and  by  Sacco  et  al.  (1990) 
using  the  FIGARO  II  balloon-borne  experiment.  However,  with  the  launch  of  the  Compton 
Gamma  Ray  Observatory,  a  positive  detection  in  the  1-30  MeV  range  has  been  reported  by  the 
COMPTEL  instrument  team  (Bennett  et  al.  1994;  Schonfelder  et  al.  1994). 

Geminga,  long  known  as  a  bright,  steady  high-energy  gamma  ray  source,  was  first  discovered 
to  be  a  0.237  second  X-ray  pulsar  by  Halpern  &  Holt  (1992)  using  ROSAT.  Once  a  period  was 
known,  pulsations  at  that  period  were  also  identified  in  EGRET  (Bertsch  et  al.  1992),  COS-B 
(Bignami  &  Caraveo  1992)  and  SAS-2  (Mattox  et  al.  1992)  gamma-ray  data.  Geminga  is  unique 
among  rotation-powered  X-ray  and  gamma-ray  pulsars  in  that,  to  date,  no  pulsed  radio  emission 
has  been  detected. 
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The  Geminga  pulse  profile  at  gamma-ray  energies,  observed  by  EGRET,  exhibits  a 
characteristic  double  peak  structure  with  peak  separation  of  ~  0.5  in  phase  (Mayer-Hasselwander 
et  al.  1994),  compared  to  a  separation  of  0.42  for  Vela  (Kanbach  et  al.  1994).  We  define  the  first 
gamma-ray  peak  to  be  the  one  preceeding  the  interpeak  emission.  Peak  FWHMs  are  a  factor  of 
2  wider  in  phase  for  Geminga  than  for  Vela. 

Measurements  of  Geminga  at  X-ray  energies  by  ROSAT  can  be  divided  into  soft  and  hard 
spectral  components,  as  described  by  Halpern  &  Ruderman  (1993).  The  lightcurve  is  energy 
dependent  with  the  single  broad  peak  observed  below  ~  0.5  keV  trailing  the  somewhat  narrower 
single  peak  observed  above  ~  0.5  keV  by  approximately  105°  in  phase. 

Evidence  exists  that  some,  if  not  most,  of  the  known  gamma-ray  pulsars  have  spectral  breaks 
somewhere  above  100  keV  (see,  for  example,  Nel  &  De  Jager  1993).  Since  the  location  of  these 
breaks  bears  on  models  of  the  emission  processes,  studies  of  Vela  and  Geminga  with  OSSE  may 
help  better  define  the  pulsar  emission  process. 

The  OSSE  results  we  report  in  this  paper  help  to  elucidate  spectral  behavior  of  both  these 
sources  via  a  detection  in  the  case  of  Vela  and  upper  limits  in  the  case  of  Geminga.  In  Section  2, 
we  discuss  the  nature  of  the  observations  made  with  OSSE.  Section  3  elaborates  on  the  analysis 
techniques  used  on  the  data  from  these  observations.  The  results  are  presented  in  Section  4  and 
discussed  in  Section  5. 


2.  Observations 


The  OSSE  instrument,  described  in  detail  by  Johnson  et  al.  (1993),  consists  of  four 
independent  phoswich  scintillation  detectors  that  operate  in  the  0.05-10  MeV  range.  Observations 
requiring  high  time  resolution  (better  than  ~  8  seconds)  are  telemetry  bandwidth  limited;  hence, 
for  both  the  Vela  and  Geminga  observations,  high  time  resolution  data  are  only  available  in 
relatively  broad  energy  bands. 

OSSE  normally  operates  in  an  on-source  off-source  chopping  mode  in  order  to  measure 
background.  However,  when  observing  a  periodically  pulsating  source,  the  flux  in  the  pulsations 
can  often  be  estimated  by  comparing  the  on-pulse  phase  regions  of  the  epoch-folded  light  curve 
to  those  off  the  pulse.  Typically,  the  pulse  profile  is  modeled  in  some  fashion.  The  phase  region 
not  included  in  the  model  pulse  is  assumed  to  be  background,  and  that  flux  level  is  subtracted 
from  the  remainder.  In  some  cases,  the  phase  region  to  be  used  for  background  determination  is 
based  on  other  criteria,  for  example,  some  previous  measurement,  perhaps  at  a  different  energy.  In 
either  case,  only  the  spectrum  of  the  pulsed  component  can  be  determined  by  OSSE  in  this  mode. 
To  measure  unpulsed  emission  from  the  source,  on-source  off-source  chopping  must  be  performed, 
reducing  on-source  livetime  substantially  (typically  by  a  factor  of  2).  The  Vela  observations 
reported  here  were  all  performed  in  “staring”  mode  in  which  no  chopping  was  performed.  Of  the 
Geminga  observations,  two  had  at  least  some  detectors  chopping,  while  one  was  performed  entirely 
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in  the  staring  mode. 

Table  1:  Vela  Pulsar  Observation  Summary 


Viewing 

Period 

Start 

Stop 

Pulsar  Livetime 
(Det-sec) 

Bands'3 

8 

8/22/91 

9/5/91 

8.08  X  105 

1,2, 3, 4 

26 / 28a 

4/23/92 

5/14/92 

5.59  X  105 

1,2, 3, 4 

301 

8/17/93 

8/24/93 

5.98  X  105 

1,2,4 

Total 

1.96  X  106 

Viewing 

Period 

Epoch  To 
JD-2440000.5 

SSB  Radio  Phase 

at  To  c 

Frequency  v 

s-1 

V 

10  “ns-2 

V 

10 -21s“3 

8 

8504 

0.034 

11.1987003711301 

-1.56791 

4.17 

26/28 

8718 

0.064 

11.1984109503280 

-1.56413 

1.49 

301 

9228 

0.403 

11.1977224439102 

-1.56079 

-0.116 

Time  resolution:  4  msec  for  all  viewing  periods 


“Vela  not  observed  from  29  April  through  6  May 
6 Band  1:  0.07  -  0.19  MeV 
Band  2:  0.22  -  0.61  MeV 
Band  3:  0.76  -  2.0  MeV 
Band  4:  2.0  -  9.7  MeV 

“Phase  of  photon  arriving  at  SSB  at  time  To  where  radio  peak  is  at  phase  0 

The  OSSE  detectors  were  used  to  observe  the  Vela  Pulsar  on  three  occasions  during  the 
two-year  interval  from  August  1991  through  August  1993.  The  individual  observations,  referred 
to  as  viewing  periods  (VP),  are  described  in  Table  1.  Note  that  the  observation  referred  to  as  VP 
26/28  was  analyzed  as  a  single  observation  even  though  it  contained  an  8-day  gap  during  which 
Vela  was  not  observed.  Live  times  are  given  in  detector-seconds  (i.e.  number  of  OSSE  detectors 
observing  times  live  time  per  detector).  The  observation  parameters  and  pulsar  ephemerides  are 
listed  in  the  table.  In  general,  band  boundaries  were  picked  in  order  to  avoid  major  background 
lines.  Note  that  the  0.77  -  2.0  MeV  band,  a  region  of  high  instrumental  background,  was  not 
included  in  VP  301.  Instead,  the  2.0  -  9.7  MeV  band  was  split  into  two  bands.  The  results  of 
these  two  bands  have  been  combined  here,  both  for  compatibility  with  the  previous  observations 
and  because  no  detection  was  made  in  the  individual  bands. 
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Table  2:  Geminga  Observation  Summary 


Viewing 

Period 

Start 

Stop 

Pulsar  Livetime  Chopped  Src  Livetime 
(Det-sec)  (Det-sec)a 

34 

7/16/92 

8/6/92 

1.22  X  106 

3.91  X  105 

310 

12/1/93 

12/13/93 

2.27  X  105 

3.75  X  105 

335b 

7/12/94 

8/1/94 

7.72  x  105 

none 

Viewing 

Epoch  Tq  EGRET  Phase 

Frequency  v  v 

Period 

JD-2440000.5 

at  Tq  c 

s'1  10-13s-2 

34,310,335  8750  0.02  4.21766909413  -1.95218  0 

Time  resolution:  8  msec  for  all  viewing  periods 


“Detector-seconds  on-source  during  chopping  observation 
^Geminga  not  observed  from  19  July  through  24  July 

“Phase  of  photon  arriving  at  SSB  at  time  To  where  EGRET  “Peak  1”  is  at  phase  0 


OSSE  observed  Geminga  on  three  occasions  from  1992  through  1994,  (see  Table  2).  As 
indicated  in  the  Table,  OSSE  operated  in  chopping  mode  during  the  VP  34  and  VP  310 
Geminga  observations,  although  two  detectors  were  staring  at  Geminga  during  VP  34.  Due  to 
the  configuration  of  the  observation,  the  chopping  mode  result  is  very  susceptible  to  systematic 
uncertainies.  We  are  continuing  evaluation  of  the  effects  and  will  report  on  the  results  of  the 
chopping  mode  observations  in  a  future  paper.  We  used  both  the  staring  mode  data  and  the 
chopping  mode  data  during  times  when  the  detectors  were  pointed  at  the  source  for  the  Geminga 
pulsar  analysis. 

During  the  Geminga  observations,  OSSE  acquired  8-ms  rate  samples  in  seven  energy  bands 
between  0.078  and  9.9  MeV.  In  order  to  improve  statistics  for  the  Geminga  pulsar  analysis,  we 
have  summed  all  the  available  data  into  two  broad  bands,  0.078  -  0.56  MeV  and  2.0  -  9.9  MeV 
(we  had  no  data  covering  the  gap  between  0.56  MeV  and  2.0  MeV). 
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3.  Analysis  and  Results 


3.1.  Epoch  Folding  Techniques 


The  pulsed  flux  analysis  consisted  of  phase- coherent  summing  (i.e.  epoch  folding)  of  4-ms 
(Vela)  or  8-ms  (Geminga)  rate  samples  to  produce  lightcurves  for  each  viewing  period  and  energy 
band.  Mean  arrival  times  at  the  solar  system  barycenter  were  computed  for  each  rate  sample 
(using  the  JPL  DE200  ephemeris)  and  the  phase  of  that  sample  was  determined  using  a  precise 
ephemeris.  In  the  case  of  Vela,  Taylor,  Nice  &  Azourmainian  (1992)  supplied  a  radio  ephemeris 
as  part  of  the  CGRO  pulsar  monitoring  program  (see  Table  1).  For  Geminga,  the  EGRET  team 
(Mattox  1994)  supplied  a  gamma-ray  based  ephemeris  which  proved  valid  over  a  wide  range  of 
epochs  (see  Table  2).  The  ephemerides  used  were  constructed  to  be  valid  for  each  of  the  viewing 
periods  involved.  Corrections  were  applied  to  account  for  a  2.042-second  clock  offset  in  the  times 
supplied  by  the  CGRO  spacecraft  during  the  first  two  Vela  observations.  Given  these  corrections, 
the  ephemerides  adequately  predicted  both  period  and  peak  absolute  phase.  We  estimate  period 
and  event  timing  uncertainties  to  be  ~  0.2  ms  for  the  determination  of  relative  phase  of  the 
beginning  of  a  given  rate  sample  and  ~  0.4  ms  for  the  absolute  phase  of  the  beginning  of  a  rate 
sample  relative  to  the  radio  or  EGRET  peak. 

We  epoch-folded  each  Vela  and  Geminga  viewing  period  separately.  For  each  rate  sample,  we 
computed  the  phase  using  the  barycentric  arrival  time  and  the  pulsar  ephemeris,  then  summed 
rate  samples  into  22-bin  phase  histograms  for  Vela  and  32-bin  phase  histograms  for  Geminga. 
Each  phase  bin  thus  represents  approximately  one  time  resolution  unit.  Since  we  had  precise 
ephemeris  information,  we  did  not  attempt  to  optimize  any  result  by  searching  frequency  space, 
nor  did  we  attempt  any  optimization  by  varying  phase  bin  widths  or  phases. 

The  epoch-folding  process  calculated  phases  relative  to  the  center  of  the  leading  radio  (Vela) 
or  gamma-ray  (Geminga)  peak  based  on  the  previously  cited  ephemerides.  For  Vela,  phase  0.0 
(the  lower  edge  of  the  first  phase  bin)  for  each  lightcurve  is  the  centroid  of  the  radio  peak.  For 
Geminga,  phase  0  is  the  nominal  centroid  of  EGRET  peak  1. 

The  analysis  produced  phase  histograms  of  counts  and  livetime  for  each  energy  band.  Since 
we  constructed  the  phase  histograms  such  that  the  absolute  phase  of  each  bin  relative  to  the 
above-mentioned  features  remained  the  same  for  all  observations,  we  were  able  to  coherently 
sum  the  histograms  from  the  three  viewing  periods.  To  characterize  lightcurve  features  with 
optimum  sensitivity,  we  summed  phase  histograms  from  several  energy  bands  together.  These 
were  normalized  by  livetime  (phase  bin  by  phase  bin)  and  energy  bandwidth  for  modeling  and 
display  purposes. 
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3.2.  Vela  Lightcurve 


We  have  detected  pulsed  emission  from  Vela  during  all  three  observations  (albeit  at  low 
statistical  significance  in  VP  301;  see  Table  3).  Figure  1  shows  the  lightcurve  summed  over  all 
the  viewing  periods  and  the  two  lowest  energy  bands  (0.07  -  0.6  MeV  overall)  .  The  horizontal 
line  represents  the  average  flux  from  the  off-pulse  or  background  phase  regions  as  described  in  the 
previous  section.  Figure  2  shows  the  OSSE  lightcurve  compared  to  other  measurements  at  lower 
and  higher  energies. 


0.0  0.2  0.4  0.6  0.8  1.0 

Phase 


Fig.  1. —  Vela  Pulsar  lightcurve  for  the  sum  of  viewing  periods  8,  26/28  and  301,  summed  over 
the  two  lowest  energy  bands  (0.07  -  0.6  MeV).  The  horizontal  solid  line  is  the  background  level 
computed  from  the  off-peak  regions.  The  dashed  line  is  the  best  fit  circular  normal  function  model 
of  the  pulse  profile. 

In  order  to  determine  on-pulse  and  off-pulse  lightcurve  regions,  we  fit  the  sum  of  the  two 
low  energy  band  lightcurves  with  a  model  consisting  of  a  constant  background  and  two  peaks 


Table  3.  Vela  Pulsar  Pulsed  Flux  Summary 


Banda 

Viewing  Period 

Peak  1  Fluxc 

Peak  2  Fluxc 

Peak  1±2  Flux' 

0.07  -  0.6b 

VP  8 

120  ±  30 

180  ±  60 

310  ±  80 

VP  26/28 

90  ±  30 

60±70 

150  ±  90 

VP  301 

50  ±40 

80  ±  70 

120  ±  90 

VP  8  +  26/28  +  301 

90  ±  20 

120  ±  40 

210  ±  50 

0.07  -  0.19 

VP  8 

190  ±  80 

240  ±  160 

430  ±  200 

VP  26/28 

180  ±  90 

60  ±  180 

240  ±  230 

VP  301 

210  ±  90 

380  ±  190 

580  ±  230 

VP  8  +  26/28  +  301 

190  ±  50 

230  ±  100 

420  ±  130 

0.22  -  0.6 

VP  8 

99  ±  31 

167  ±  65 

267  ±  81 

VP  26/28 

61  ±  37 

62  ±  76 

123  ±  95 

VP  301 

-7±  37 

-19  ±  76 

-26  ±  95 

VP  8  +  26/28  +  301 

57  ±  20 

83  ±41 

140  ±  52 

0.76  -  2.0 

VP  8 

-8  ±  14 

-30  ±  30 

-38  ±  -38 

VP  26/28 

7±  16 

-6  ±  34 

0.7  ±  43 

VP  301 

VP  8  +  26/28 

-2  ±  11 

-20  ±  22 

-22  ±  28 

2.0  -  9.7 

VP  8 

0.7  ±  1.5 

5.6  ±  3.1 

6.2  ±  3.9 

VP  26/28 

1.7  ±  1.8 

0.6  ±  3.7 

2.3  ±4.6 

VP  301 

-1.4  ±  1.8 

5.1  ±  3.6 

3.6  ±4.5 

VP  8  +  26/28  +  301 

0.3  ±  1.0 

4.0  ±  2.0 

4.3  ±  2.5 

aBands  in  MeV 

bNote  that  this  is  the  concatenation  of  the  next  two  bands 
c10-6  photons/cm2-s-MeV,  phase  averaged 
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Phase 


Fig.  2. —  Vela  Pulsar  lightcurves  from  ROSAT  (Ogelmanet  al.  1993),  OSSE,  COMPTEL  (Bennett 
et  al.  1994)  and  EGRET  (Kanbach  et  al.  1994),  aligned  in  phase  with  the  radio  phase  at  0.0.  The 
dotted  lines  show  the  EGRET  peak  regions. 

represented  by  circular  normal  functions.  The  circular  normal  function  is  defined  by  the  relation: 

I p  ^  k  cos(2w(x  —  xr.)) 

~  Io(k) 

where  A  controls  the  peak  amplitude,  k  is  the  “compactness”,  which  is  inversely  related  to  the 
peak  width,  xc  is  the  phase  of  the  center  of  the  peak,  and  Iq  is  the  modified  Bessel  function  of 
order  zero.  The  circular  normal  function  has  the  property  that,  although  it  is  roughly  gaussian  in 
shape,  it  is  periodic  in  phase  with  period  1.  Hence,  it  is  well  suited  to  modeling  lightcurve  peaks 
that  “wrap  around”  the  phase  histogram.  All  parameters  in  the  peaks  and  background  were  free 
to  vary.  We  used  these  models  to  determine  peak  positions  and  widths,  but  not  peak  amplitudes, 
as  described  below.  The  resulting  peak  positions  and  widths  are  listed  in  Table  4. 

We  defined  the  “peak”  or  “on-source”  phase  regions  as  the  regions  of  the  light  curve  within 
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Table  4:  Vela  Pulsar  Pulse  Profile  Model  Best  Fit  Parameters3 


Peak  1  Phase  O.ll+ooi 

Peak  Separation  0.44  ±  0.06 

Peak  1  FWHM  0.08to!o4 

Peak  2  FWHM  0.22±g;?| 


“All  parameters  are  in  units  of  folded  lightcurve  phase  (0.0  -  1.0) 

the  full  width  at  10%  maximum  of  the  two  best-fit  peak  models.  Since  we  saw  no  evidence  for 
any  other  significant  excesses,  the  remainder  of  the  lightcurve  was  treated  as  unpulsed  emission 
and  used  to  determine  background  for  spectral  analyses.  For  each  energy  band,  the  background 
regions  in  the  original  count  and  livetime  histograms  were  summed,  normalized  by  the  ratio  of  the 
livetime  from  the  peak  region  to  that  from  the  background  region,  and  subtracted  from  the  total 
counts  from  the  peak  region.  The  resulting  difference,  still  in  counts,  was  then  normalized  by  the 
total  live  time  for  the  entire  lightcurve,  resulting  in  a  phase-averaged  peak  rate.  This  technique  is 
preferable  to  using  the  circular  normal  function  peak  amplitudes  to  provide  peak  fluxes,  since  it  is 
not  as  vulnerable  to  misrepresentation  of  the  data  by  the  model  and  is  more  straightforward  for 
uncertainty  calculation. 

The  significance  of  the  detection  was  determined  using  the  maximum  likelihood  ratio  technique 
described  by  Li  &  Ma  (1983)  applied  to  the  summed  fluxes  in  the  “source”  and  “background” 
regions.  The  source  region  is  defined  to  be  the  sum  of  the  two  peak  regions,  while  the  background 
region  is  everything  else.  The  distribution  of  counts  in  the  background  region  can  be  represented 
by  expectation  value  {Ng)  while  the  distribution  of  counts  in  the  source  region  is  {Ns)  +  a{Ng), 
where  {Ns)  is  the  expectation  value  of  the  number  of  excess  pulsed  counts  from  the  source 
and  a  =  ts/tg  is  the  ratio  of  source  region  to  background  region  livetimes.  The  test  computes 
a  ratio  of  the  maximum  likelihoods  in  which  the  numerator  is  the  likelihood  of  obtaining  the 
observed  results  given  the  null  hypothesis  {Ns)  =  0.  The  denominator  is  the  maximum  likelihood 
of  obtaining  the  observed  results  given  that  {Ns)  can  be  non-zero.  Li  &  Ma  have  shown  that 
S  =  v  —  2  In  A  (where  A  is  the  likelihood  ratio)  is  distributed  as  the  absolute  value  of  a  standard 
normal  variable  (i.e.  zero  mean  and  unit  variance).  The  probability  (as  shown  in  Table  5)  of 
obtaining  a  value  of  S  equal  to  or  greater  than  the  one  observed  tests  the  null  hypothesis.  A  small 
probability  means  that  the  null  hypothesis  can  be  rejected  with  some  confidence,  i.e.  that  there  is 
a  significant  excess  (or  deficiency)  of  counts  in  the  pulse  region. 

Table  5  displays  the  matching  results  of  the  Maximum  Likelihood  Ratio  (MLR)  test  for  each 
peak  and  viewing  period  summed  over  the  two  lowest  energy  bands. 

The  overall  MLR  significance  of  detection  of  the  first  gamma  ray  pulse  profile  peak  integrated 
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Table  5:  Vela  Pulsar  Maximum  Likelihood  Ratio  (MLR)  Test  Summary 

Viewing  Period  Peak  1  MLR  Peak  2  MLR  Peak  1+2  MLR 


VP  8 

2.9  X  10“5 

1.3  x  10-3 

VP  26/28 

6.2  x  10-3 

2.0  x  10-1 

VP  301 

1.3  X  10-1 

1.7  x  10-1 

VP  8  +  VP  26/28  +  VP  301 

2.4  x  10“6 

1.7  x  10-3 

4.0  x  10“5 
4.9  x  10-2 

1.1  x  10-1 

2.1  x  10“5 


over  all  three  observations  is  ~  4.6<r  in  the  sum  of  the  two  lower  energy  bands.  We  detect  the 
second  gamma  ray  peak  at  only  marginal  significance  (~  3<r)  in  the  same  band.  Although  Table  5 
gives  the  impression  that  we  did  not  detect  the  pulsar  at  all  during  VP  301,  note  that  this  table 
represents  data  in  the  sum  of  the  two  lower  energy  bands  (0.07  -  0.6  MeV).  Comparing  detections 
in  the  lowest  band  alone  (0.07  -  0.19  MeV)  in  Table  3,  we  see  that  VP  301  did  indeed  contribute 
to  the  overall  detection  significance. 

We  see  no  significant  evidence  of  temporal  variability  in  either  peak  or  any  of  the  energy 
bands.  Table  6  shows  the  results  of  computing  x2  assuming  a  constant  flux  model  (with  flux 
equal  to  the  mean)  for  each  peak.  For  none  of  the  combinations  of  peaks  and  bands  can  the 
constant  flux  (i.e.  no  temporal  variability)  model  be  rejected  with  as  much  as  95%  confidence  (the 
probability  of  observing  the  data  with  a  constant  model  assuming  random  fluctuations  is  given  in 
parentheses  in  the  table).  Although  the  source  is  not  detected  in  the  0.22  -  0.6  MeV  band  during 
VP  301  (note  the  higher  x2  values  in  that  band),  the  resulting  flux  decrease  is  not  significant, 
as  indicated  by  the  probabilities  of  random  fluctuations  producing  the  observed  x2 .  In  addition, 
we  have  compared  lightcurves  for  each  band  between  pairs  of  viewing  periods  (subtracting  the 
background  in  each  case)  and  see  no  evidence  for  significant  variability.  For  the  remainder  of  this 
discussion,  we  will  treat  the  sum  of  the  three  viewing  periods  only. 

The  lightcurve  looks  qualitatively  similar  to  the  lightcurves  measured  at  higher  energies 
(Buccheri  et  al.  1978;  Grenier  et  al.  1988;  Kanbach  et  al.  1994).  We  have  characterized  our 
observed  lightcurve  shape  via  a  model  as  mentioned  above.  The  results  of  this  characterization  (for 
the  optimal-sensitivity  0.07  -  0.60  MeV  band)  are  displayed  in  Table  4  and  plotted  in  Figure  1. 
OSSE  measures  a  first  peak  width  of  O.O8+Q4  in  phase  in  the  0.07  -  0.60  MeV  band.  The  EGRET 
first  peak  is  ~  0.03  wide  in  phase,  which,  although  narrower,  is  not  statistically  inconsistent  with 
the  OSSE  result.  The  second  peak  appears  to  be  almost  a  factor  of  three  broader  than  the  first, 
as  opposed  to  EGRET  observations,  which  indicate  that  the  second  peak  is  ~  1.7  times  broader 
than  the  first  above  100  MeV  (Kanbach  et  al.  1994).  However,  not  only  is  the  width  of  the  second 
peak  not  well  determined  (as  indicated  in  Table  4),  but  the  F-test  for  models  with  and  without 
the  second  peak  indicates  that  it  is  only  required  at  the  ~  2.6<r  level.  Hence,  although  we  have 
an  indication  that  the  second  peak  is  broader  than  the  first,  we  cannot  make  this  claim  with 
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Table  6:  Vela  Pulsar  Pulsed  Flux  Variability  Test  Summary 

Banda  x2/dof  Peak  lb  x2/dof  Peak  2  %2/dof  Peak  1+2 


0.07  -  0.6C 
0.07  -  0.19 
0.22  -  0.6 
2.0  -  9.7 


1.51  (0.22) 
0.02  (0.98) 
2.41  (0.09) 
0.82  (0.44) 


1.15  (0.32) 
0.77  (0.47) 
1.84  (0.16) 
0.63  (0.54) 


1.61  (0.20) 
0.56  (0.57) 
2.78  (0.06) 
0.23  (0.79) 


“Bands  in  MeV 

6Tests  hypothesis  that  individual  viewing  periods  are  randomly  distributed  about  the  mean  of  all  the  viewing 
periods;  number  in  parentheses  is  Probability(%2  >  xlbsvd) 

“Note  that  this  is  the  concatenation  of  the  next  two  bands 


confidence.  The  positions  of  the  OSSE  peaks,  as  represented  by  the  best  fit  model,  are  consistent 
with  the  EGRET  peak  postions  given  in  Kanbach  et  al.  (1994). 


3.3.  Geminga  Lightcurve 


We  have  analyzed  the  Geminga  data  for  pulsed  emission  in  a  fashion  similar  to  Vela,  using 
ephemerides  provided  by  the  EGRET  team  (Mattox  1994).  We  saw  no  evidence  of  pulsations  in 
either  the  low  or  high  energy  bands  in  the  individual  observation  lightcurves  nor  in  the  lightcurve 
representing  the  sum  of  the  three  observations.  The  latter  are  shown  in  Figure  3.  The  EGRET 
30  -  100  MeV  lightcurve,  also  shown  in  this  figure,  is  taken  from  Mayer-Hasselwander  et  al. 
(1994).  It  has  been  phase  aligned  with  the  OSSE  lightcurve  using  the  relative  epochs  and  phases. 
Note  that  the  phase  axis  has  been  shifted  such  that  the  zero  of  phase,  which  in  this  representation 
corresponds  to  the  centroid  of  the  first  EGRET  peak2,  is  in  the  middle  of  the  figure. 

We  can  demonstrate  in  several  ways  that  there  is  no  statistically  significant  evidence  for 
peaks  in  the  Geminga  lightcurve.  First,  if  we  simply  compute  the  pulsed  flux  in  the  EGRET 
peak  regions,  using  for  background  the  phase  regions  not  included  in  the  EGRET  peaks,  we  see 
no  evidence  of  significant  excess.  Further,  assuming  that  the  data  are  represented  by  the  mean, 
we  calculate  %2.  For  the  low  energy  band,  %2/dof=  1.33  for  31  dof.  The  probability  of  observing 
this  value  or  higher  given  only  random  fluctuations  about  the  constant  mean  is  0.10.  A  similar 
calculation  for  the  high  band  gives  %2/dof=  0.91  with  a  probability  of  0.60.  Hence  we  cannot 


2The  phase  axis  is  referenced  to  the  nominal  EGRET  peak  as  defined  by  the  EGRET  ephemeris.  This  may  differ 
slightly  from  the  peak  displayed  in  the  figure,  since  the  peak  position  varies  slightly  with  energy 
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Phase 


Fig.  3. —  Geminga  lightcurves  for  the  sum  of  VP  34,  VP  310  and  VP  335  and  two  energy  bands, 
together  with  30-100  MeV  EGRET  lightcurve  from  Mayer-Hasselwander  et  al.  (1994). 

reject  the  hypotheses  that  the  parent  distributions  are  independent  of  phase. 

Next,  we  have  computed  the  Li  &  Ma  (1983)  MLR  for  the  positions  of  each  EGRET  peak.  In 
no  case  (i.e.  peak  1,  peak  2,  or  both  peaks  in  either  energy  band)  can  we  reject  the  null  hypothesis 
of  zero  source  counts  with  better  than  1.5<r  (87%)  confidence.  In  addition,  we  have  examined  the 
first  half  of  the  EGRET  interpulse  2  (12)  region  (the  interpulse  region  with  less  flux),  since  this 
region  has  the  softest  EGRET  spectrum  and  the  first  half  of  the  region  appears  to  be  the  source  of 
the  flux  in  the  30-100  MeV  light  curve.  We  cannot  reject  the  null  hypothesis  for  this  phase  region 
to  any  better  than  0.6<r. 

ROSAT  observations  of  pulsed  emission  from  Geminga  (Halpern  &  Holt  1992;  Halpern  & 
Ruderman  1993)  have  demonstrated,  in  the  0.5  -  1.5  keV  band,  a  peak  in  the  lightcurve  between 
EGRET  peaks  1  and  2.  Halpern  (1995)  has  indicated  that  the  absolute  phase  of  the  ROSAT 
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lightcurve  relative  to  the  EGRET  lightcurve,  as  shown  in  Figure  2  of  Halpern  &  Ruderman  (1993) 
is  not  reliable.  Therefore,  we  have  not  attempted  to  phase  align  this  result  with  the  OSSE  light 
curves. 


3.4.  Vela  Spectrum 


We  used  a  forward-folding  process  to  compare  the  data  to  various  model  spectra.  In  this 
procedure,  a  model  spectrum  was  folded  through  the  instrument  response  and  the  resulting  model 
count  rate  spectrum  was  compared  to  the  data  in  a  x2  sense.  The  model  spectra  and  responses 
were  created  and  multiplied  in  narrow  energy  bands,  then  integrated  over  the  broad  OSSE  bands 
so  that  spectral  variations  over  each  band,  appropriately  weighted  by  response,  were  accounted  for. 
The  fitting  routine  used  was  based  on  the  Marquardt  algorithm  with  data  points  weighted  by  their 
uncertainties.  The  routine  used  no  special  treatment  for  upper  limits  (i.e.  insignificant  data  points 
were  treated  like  all  the  rest).  The  model  parameters  were  varied  until  a  model  that  generates  a 
minimum  x2  was  determined.  Since  our  background  counting  rates  were  high  and  energy  bands 
broad,  this  was  also  guaranteed  to  be  a  maximum  likelihood  estimate.  Parameter  uncertainties 
were  determined  by  2 -mapping  using  a  procedure  based  on  that  described  by  Lampton,  Margon 
&  Bowyer  (1976). 

Data  from  other  experiments  (e.g.  EGRET  and  COMPTEL)  were  fitted  simultaneously 
with  OSSE  spectra  using  the  same  forward  folding  technique.  However,  since  we  did  not  possess 
instrument  responses  for  these  data  sets,  published  photon  spectra  were  used  and  unit  responses 
assumed.  The  model  spectrum  was  still  integrated  over  each  band,  but  weighting  by  response 
within  a  band  was  not  possible,  so  the  results  are  somewhat  approximate  and  should  be  viewed 
qualitatively. 

We  generated  photon  spectra  from  the  data  by  multiplying  each  data  point  in  the  count 
spectrum  by  the  ratio  of  the  best-fit  model  photon  spectrum  to  the  best-fit  model  count  spectrum. 
The  resulting  spectrum  was  dependent  on  the  model  chosen,  although  in  the  case  of  the  Vela 
Pulsar,  the  uncertainty  in  the  result  due  to  model  dependence  was  small  relative  to  the  statistical 
uncertainty.  The  resulting  photon  spectra  are  compiled  in  Table  3. 
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10'4  10~2  10°  102  104 


Energy  (MeV) 

Fig.  4. —  Vela  Pulsar  spectrum  for  the  sum  of  viewing  periods  8,  26/28  and  301.  Flux  is  the  sum 
of  both  peaks,  background  subtracted  and  averaged  over  the  entire  lightcurve.  It  is  presented  as 
(E)xphoton  flux.  Upper  limits  are  2<r  above  zero.  The  data  and  models  displayed  are:  EGRET 
(Kanbach  et  al.  1994);  COMPTEL  (Schofelder  et  al.  1994);  UCR  (Turner  et  al.  (1984);  COS-B 
1981  best  fit  two  power  law  model  (Grenier  et  al.  1988);  ROSAT  (Ogelman  et  al.  1993);  Vela  Outer 
Gap  Model  (Cheng  et  al.  1986);  and  Polar  Cap  Models  (DH,  Daugherty  &  Harding  1995;  and  SDM, 
Sturner  et  al.  1995).  The  Vela  Outer  Gap  Model  has  been  adjusted  in  shape  and  normalization  to 
fit  the  data,  while  the  Polar  Cap  Models  have  been  normalized  only.  The  shaded  region  represents 
the  actual  COS-B  energy  range,  while  the  dashed  lines  extrapolate  the  COS-B  spectrum  to  lower 
energies. 
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Figure  4  displays  the  sum  of  the  first  and  second  peak  emission  (i.e.  the  difference  of  peak 
region  and  off-peak  or  background  region  rates).  In  order  to  compress  the  display  in  the  OSSE 
energy  range,  the  spectrum  is  plotted  as  (E)xphoton  flux  (i.e.  an  energy  flux),  where  (E)  is  defined 
as  the  mean  band  energy,  computed  by  averaging  the  energy  weighted  by  a  representative  spectrum 
model  over  each  band.  Note  that  we  use  this  scheme  strictly  for  dynamic  range  compression  in 
the  plot.  All  of  our  models  are  specified  as  photon  spectra,  according  to  common  practice  in  this 
energy  range.  Fluxes  are  normalized  by  division  by  the  total  livetime  and  are  hence  time  averaged 
rather  than  instantaneous  values.  Data  points  less  than  one  standard  deviation  above  zero  have 
been  plotted  as  2<r  upper  limits  (that  is,  two  standard  deviations  above  zero). 


(10  4  Photons/cm2-s-MeV  at  100  keV) 

Fig.  5. —  Confidence  contours  for  a  single  power-law  fit  to  the  OSSE  Vela  pulsar  spectrum  (peaks 
1  +  2).  Confidence  intervals  and  number  of  interesting  parameters  are  indicated  on  the  contours. 
Intervals  are  determined  using  the  prescription  of  Lampton,  Margon  &  Bowyer  (1976) 

The  spectrum  is  rather  hard,  particularly  below  0.60  MeV.  The  OSSE  spectrum  is  well- 
represented  by  a  single  power  law  (see  Figure  5  for  parameter  confidence  contours).  The  best  fit 
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power  law  index  for  the  sum  of  both  peaks  is  —1.3  ±  0.2,  while  the  best  fit  index  for  peak  1  is 
—  1.5^0  3  and  for  peak  2  is  —1.1  ±  0.2.  Note  that  the  uncertainty  for  peak  2  is  comparable  to  that 
for  peak  1  even  though  peak  1  appears  more  significant  in  Table  5.  This  is  because  the  Table  is 
computed  for  the  sum  of  the  two  lower  energy  bands,  where  peak  1  is  strongest,  while  a  significant 
portion  of  the  peak  2  flux  comes  from  the  highest  energy  band  (see  Table  3),  which  is  not  included 
in  the  computation  of  MLR.  The  peak  1  and  peak  2  results  are  consistent  with  each  other  to 
within  statistics.  The  OSSE  result  thus  contains  no  evidence  of  phase- dependent  spectral  shape. 
We  will  discuss  several  other  models  for  the  data  in  Section  4. 


3.5.  Geminga  Spectrum 


In  the  case  of  Geminga,  in  which  no  significant  pulsations  were  observed,  we  used  the  following 
formula  from  Ulmer  et  al.  (1991),  to  compute  upper  limits: 

NVC^t  I  P 

A E  Aeff  t  V  1  -  p 

where  N  is  the  confidence  level  for  the  upper  limit  in  units  of  standard  deviation  (taken  to  be  2), 
P  is  the  pulsar  duty  cycle,  Ct„t  is  the  total  number  of  counts  in  the  lightcurve,  A E  is  the  energy 
band,  Aeff  the  effective  area  integrated  across  that  band  and  t  the  livetime.  For  a  pulsar  duty 
cycle,  we  have  used  0.5  (the  most  conservative  choice).  Although  the  EGRET  peaks  are  much 
narrower  than  this,  ROSAT  Halpern  &  Ruderman  (1993)  and  COMPTEL  (Kuiper  et  al.  1995) 
report  different  light  curve  shapes  and  peak  locations,  making  prediction  of  the  light  curve  shape 
in  the  OSSE  band  difficult. 

Figure  6  displays  the  OSSE  upper  limits  computed  as  described  above,  together  with  a 
spectrum  measured  by  COMPTEL  (Kuiper  et  al.  1995)  and  the  EGRET  total  pulsed  spectrum 
together  with  an  extrapolation  of  the  EGRET  best  fit  power  law  model  to  lower  energies 
(Mayer-Hasselwander  et  al.  1994). 

Although  not  plotted  in  the  figure,  we  have  also  considered  the  OSSE  upper  limits  for  the 
first  half  of  the  EGRET  12  region  as  discussed  above.  These  are  lower  than  the  plotted  limits  by 
the  ratio  of  the 

fb 

terms  in  the  upper  limit  expression,  where  this  term  is  equal  to  1  for  /3  =  0.5  and  is  equal  to  0.5 
for  P  =  0.2,  the  duty  cycle  for  the  first  half  of  the  12  region.  Hence  the  upper  limits  for  this  phase 
assumption  are  lower  by  a  factor  of  2  than  those  shown. 
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Fig.  6. —  Geminga  spectral  upper  limits  from  OSSE  and  COMPTEL,  together  with  an 
extrapolation  of  the  EGRET  best  fit  power  law  model.  Also,  the  reported  pulsed  flux  from 
FIGARO  II  (Sacco  et  al.  1994). 

4.  D  iscussion 

The  Vela  lightcurve  as  measured  by  OSSE  is  roughly  similar  to  the  lightcurve  observed  at 
higher  energies  (see  Figure  2).  However,  we  observe  weak  evidence  that  the  second  peak  could 
be  somewhat  broader  than  that  observed,  for  example,  by  EGRET.  Hence,  the  OSSE  data  may 
serve  as  a  bridge  between  the  high  energy  data  and  the  lightcurve  measured  at  X-ray  energies  by 
ROSAT  (Ogelman  et  al.  1993).  The  latter  has  a  broad,  complex  pulse  centered  at  a  phase  of 
0.7  relative  to  the  radio  pulse.  A  broad  second  pulse  in  the  X-ray  band  has  been  predicted  by 
Romani  &  Yadigaroglu  (1994). 

The  emerging  picture  of  the  Vela  pulsed  spectrum  consists  of  a  hard  spectrum  at  low 
gamma-ray  energies,  breaking  to  a  rather  softer  spectrum  at  higher  energies  (Figure  4).  The 
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EGRET  and  COMPTEL  spectra  are  total  pulsed  flux  in  each  case.  Fitting  a  power  law  model  to 
OSSE,  COMPTEL  (Schonfelder  et  al.  1994)  and  EGRET  (Kanbach  et  al.  1994)  data  results  in 
an  unacceptably  high  %2/dof  (39  for  16  degrees  of  freedom).  Using,  instead,  a  broken  power  law 
model  requires  a  break  between  the  two  power  law  components  at  26  ±  8  MeV.  The  broken  power 
law  model  results  in  a  %2/dof  of  7.7  for  14  degrees  of  freedom,  better  than  the  single  power  law 
but  still  not  very  good.  Since  the  actual  spectral  shape  may  roll  over  smoothly  rather  than  with 
an  abrupt  break,  the  break  energy  of  the  broken  power  law  is  only  a  rough  indication  of  spectrum 
behavior  and  may  be  dependent  on  the  energy  range  used  in  the  fit.  If  we  perform  the  same 
fit  to  OSSE,  COMPTEL  and  EGRET  peaks  1  and  2  data  separately,  relatively  poor  OSSE  and 
COMPTEL  statistics  do  not  allow  us  to  rule  out  the  same  break  energy  and  magnitude  to  each 
peak,  or  to  the  sum  of  the  peaks. 

The  X-ray  spectrum  displayed  in  Figure  4  is  the  single  blackbody  model  that  best  represents 
the  pulsed  emission  observed  by  ROSAT  (Ogelman  et  al.  1993).  There  has  been  no  report  of  a  hard 
component  in  the  pulsed  emission,  although  such  a  component  has  been  suggested  by  Ogelman 
(1993)  for  the  total  emission  from  the  point  source.  Bridging  the  gap  between  the  ROSAT  and 
OSSE  observations  could  confirm  the  details  of  the  magnetospheric  gamma  ray  emission  at  low 
energies,  where  modeling  is  difficult.  This  should  be  a  high  priority  for  observations  of  Vela  during 
future  missions. 

Figure  4  displays  three  theoretical  model  spectra  for  pulsed  emission  from  Vela.  The  Vela 
outer  gap  model  (Cheng,  Ho  &  Ruderman  1986  and  Ho  1993),  and  two  polar  cap  emission  models 
(Daugherty  &  Harding  1995,  labeled  DH;  and  Sturner,  Dermer  &  Michel  1995,  labeled  SDM) 
represent  the  data  with  varying  degrees  of  success.  The  only  model  with  a  significant  disagreement 
in  the  OSSE  band  is  the  last,  which  rolls  over  to  a  spectrum  that  is  too  hard  below  1  MeV  and 
too  soft  in  the  EGRET  energy  range.  We  emphasize  that  the  poor  agreement  with  the  data  of  the 
two  polar  cap  models  compared  to  the  two  outer  gap  model  is  somewhat  misleading,  since  the 
outer  gap  model,  an  analytic  function  with  three  free  parameters,  was  actually  fit  to  all  the  data, 
while  the  other  two,  available  only  in  tabular  form,  were  normalized  only. 

The  Vela  outer  gap  model,  as  used  here,  has  two  parameters  in  addition  to  normalization. 
These  are  Emax,  the  maximum  secondary  synchrotron  photon  energy,  and  unun  =  7mm/7max, 
where  7 max  is  the  maximum  energy  achieved  by  electrons  in  the  gap  accelerator,  and  7 m8„  is  the 
low  energy  cutoff  in  the  electron  spectrum  caused  by  electrons  escaping  across  the  light  cylinder 
before  losing  all  their  energy  to  synchrotron  radiation.  The  radiation  spectrum  is  a  smoothly 
varying  power  law  from  Emin  =  u2nnn  x  Ernax  to  Ernax.  Below  Emin,  the  spectrum  should  harden, 
with  E-2/3  the  hardest  possible  spectrum  expected  (Ho  1993).  Fitting  the  outer  gap  model  to 
OSSE,  COMPTEL  and  EGRET  data  results  in  Emax  =  41  ±  12  GeV  and  um8„  =  0.021  ±  0.008 
which  leads  to  Ema„  =  18  MeV.  The  latter  is  consistent  with  the  break  energy  found  using  a 
broken  power  law  model.  A  simple  power  law  fit  to  the  OSSE  data  alone  has  a  best  fit  index  of 
—  1.3  ±  0.2,  which  is  softer  than  the  hardest  predicted  index  of  —2/3  and,  hence,  does  not  reject 
these  models  on  the  basis  of  exceeding  the  maximum  hardness.  Note  that  the  best  fit  Vela  outer 
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gap  model  does  not  represent  the  high  energy  cutoff  present  in  the  EGRET  spectrum  at  ~  4  GeV. 
This  is  not  surprising  given  the  best-fit  value  for  Emax  of  41  GeV.  If  Emax  is  lowered  to  produce  a 
cutoff  at  the  lower  energy,  the  modeled  cascade  produces  insufficient  photons  in  the  OSSE  energy 
range. 

Turner  et  al.  (1984)  have  reported  a  detection  of  pulsed  emission  from  Vela  in  the  range 
0.3-30  MeV  based  on  a  balloon  flight  in  November  1981.  Their  results  are  shown  in  Figure  4. 
Their  reported  fluxes  are  significantly  higher  and  perhaps  softer  than  those  measured  by  the 
CGRO  instruments;  however,  contemporaneous  measurements  by  COS-B  in  1981  (Grenier  et  al. 
1988)  were  also  soft  compared  to  the  EGRET  measurements,  which,  according  to  Kanbach  et  al. 
(1994),  have  shown  no  evidence  of  variability  from  May  1991  through  November  1992.  Figure  4 
shows  that  an  extrapolation  of  the  COS-B  spectrum  below  300  MeV  is  still  marginally  too  low 
to  explain  the  UCR  result,  although  a  downward  extrapolation  of  the  spectrum  above  300  MeV 
would  intersect  the  Turner  result  reasonably  well.  The  coincidence  of  both  1981  measurements 
being  higher  than  the  CGRO  results  hints  at  long  term  variability  in  the  Vela  pulsed  flux. 

The  Geminga  spectrum  (Figure  6)  in  low  energy  gamma  rays  is  not  as  well  determined  as  for 
Vela.  The  OSSE  2<r  upper  limit  from  0.08  -  0.6  MeV  is  low  compared  to  the  extrapolation  of  the 
EGRET  best  fit  spectrum,  as  shown  here  with  an  extrapolated  68%  confidence  region  (Mayer- 
Hasselwander  et  al.  1994).  However,  the  disagreement  is  not  very  significant.  If  we  integrate 
the  EGRET  extrapolation  of  the  lower  OSSE  band,  it  disagrees  with  the  OSSE  null  detection 
by  1.5cr,  a  result  almost  entirely  driven  by  the  uncertainty  in  the  EGRET  extrapolation.  The 
COMPTEL  result  is  also  more  or  less  consistent  with  the  EGRET  extrapolation.  Therefore,  we 
cannot  insist  on  a  break  in  the  spectrum  similar  to  that  observed  in  Vela.  However,  the  spectrum 
must  break  somewhere  in  order  not  to  exceed  the  ROSAT  flux.  The  break  cannot  be  much  below 
10  keV  without  requiring  an  unreasonably  hard  spectrum  below  the  break,  and  it  cannot  break 
much  above  3  MeV  without  contradicting  the  COMPTEL  result.  The  OSSE  upper  limit  disagrees 
strongly  with  the  claimed  detection  by  FIGARO  II  (Sacco  et  al.  1994). 


5.  Summary 


OSSE  observations  of  the  Vela  and  Geminga  Pulsars  have  detected  the  former  and  set 
significant  upper  limits  on  the  spectrum  of  the  latter  in  the  low  energy  gamma  ray  range.  In 
both  cases,  OSSE  (and  COMPTEL)  results,  when  extrapolated  up  in  energy,  require  a  break  to  a 
softer  spectrum  in  order  to  intersect  EGRET  data  above  70  MeV.  For  Vela,  this  break  appears  in 
the  20-30  MeV  range,  while  for  Geminga  it  is  not  well  determined.  The  Vela  Pulsar  light  curve 
measured  by  OSSE  is  similar  to  that  at  higher  energies,  within  the  limited  statistics  of  the  OSSE 
detection.  The  Vela  pulsed  spectrum  in  the  OSSE  range  is  quite  hard,  although  not  so  hard  that 
it  violates  constraints  of  synchrotron  production  models.  The  data  are  well  represented  by  an 
outer  gap  model,  but  polar  cap  models  are  not  firmly  excluded. 
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